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We calculate the cross section for diffractive photoproduction of opposite-charge pseudoscalar 
meson pairs M + M~ — ty + tt' ~, K+K~, D + D~ and B + B~ in a broad range of center-of-mass 
energies relevant for GlueX/Hall D, FOCUS, COMPASS and HERA experiments. In the case of 
tt + tv~ production we find that the interference of the p° resonance and the two-pion continuum 
leads to a considerable deformation of the shape of p° in agreement with the data from the ZEUS 
collaboration. We also discuss the spectral shape of the p° as a function of the momentum transfer 
and the contribution of higher partial waves to the -k + tt~ mass spectrum. We predict a sizeable 
energy-dependent forward-backward asymmetry in the Gottfried- Jackson frame. For the heavy 
meson production we find that the cross section for diffractive production increases much slower 
than the one for open charm or bottom production. We discuss lower and upper limits for the cross 
sections for diffractive production of D + D~ and B + B~ pairs, which we find can be as large as 10% 
of the open flavor production. 
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I. INTRODUCTION 

The mechanism of diffractive photoproduction of 
charged pion pairs has been proposed long ago 0, 0, 0, 
QLIMlI' At that time there, however, no good data was 
available for constraining details of the model. In recent 
years a body of data from ttN and KN scattering has 
been collected (see for instance |7()- In particular the 
ZEUS and HI collaborations at HERA have measured 
diffractive production of the p° meson. Most of the the- 
oretical effort, however, concentrated on the description 
of the total cross section in different pQCD inspired ap- 
proaches and not on the description of the two-pion con- 
tinuum invariant mass distribution. 

In real photoproduction the ZEUS collaboration ob- 
served a strong asymmetry in the two-pion invariant mass 
spectrum around the peak position [8j of the p°. In elec- 
troptoduction the asymmetry seems to decrease with in- 
creasing photon virtuality |9(, Q 2 and momentum trans- 
fer t' = \t — t m i n \ [lCj. One would expect, whatever the 
non-resonant mechanism might be, that it will also be 
present at small Q 2 and small t' and thus needed to iso- 
late the p° production cross section. 

The GlueX (Hall D) project at the Jefferson Lab will 
study mesonic resonances focusing on gluonic excitations. 
The diffractive production of charged meson continuum 
may produce large backgrounds to the channels of inter- 
est. Recently exotic J PC = 1~ + candidates have been 
reported by the E852 [13 and the Crystall Barrel [13 col- 
laborations. Candidates for the exotic states are rather 



broad |lll 1 1 121] and may indeed have a large component 
originating from production of the meson continuum pro- 
duction |lall4| . 

Recently the FOCUS collaboration at Fermilab found 
a new state in the K + K~ final state at 1.75 GeV 0] 
wich may also require a good understanding of the KK, 
continuum production. 

In the present paper we study the tt + tt~ and K + K~ 
channels. In this case model ingredients are well con- 
strained by the ttN and KN data. The invariant 
mass distribution and possible 7r + — 7r~ polar angle asym- 
metries are presently being analyzed by the COMPASS 
collaboration [l^. We also compare the results of our 
calculation with the experimental data at higher energies 
from HERA and medium energy relevant for the future 
GlueX/Hall D experiment at TJNAF. 

The inclusive production of heavy charmed mesons in 
electro- or pho-production off proton is routinely used to 
study the gluon distribution in the nucleon. The stan- 
dard QCD approach is based on the production of heavy 
quark - heavy antiquark pairs at the parton level from 
photon-gluon fusion followed by a fragmentation to heavy 
flavored hadrons. The formalism which we present for 
diffractive production of pairs of light charged mesons 
should be also valid for the production of pairs of heavy 
mesons, D + D~ or even B + B~ . This may be interesting 
in the context of a deficit in the open bb production in 
photon-proton [2(j and photon-photon [2lJ collisions. 

Recently the FOCUS collaboration has analyzed the 
azimuthal correlations between D D mesons [221. It was 
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pointed out very recently |23j ] that heavy meson correla- 
tions are very useful to study unintegrated gluon distri- 
butions in the nucleon. To the best of our knowledge, 
the contribution of the diffr active mechanism, discussed 
here, has not been estimated in this context. It is also in- 
teresting to investigate how large the diffractive produc- 
tion of B + B~ pairs might be, compared to the standard 
pQCD mechanism of bb production discussed above. In 
photoproduction, one would naively expect a relative en- 
hancement of the ratio of diffractive B + B~ (charge 1) to 
the standard bb (charge 1/3) production as compared to 
the ratio of the D + D~ (charge 1) to the cc (charge 2/3) 
production by a factor of 4 in the cross section. A better 
understanding requires, however, more detailed insight 
into the dynamics of the process. In the present paper 
we shall present estimates of such contributions. In par- 
ticular, we discuss several aspects of the opposite-charge 
pseudoscalar continuum in photoproduction. 



We note that a similar model has been recently applied 
to describe 7T7T and KK photoproduction in Ref . . 



II. MODEL OF THE CONTINUUM 

The dominant mechanisms of diffractive production of 
opposite-charge meson pairs are shown in Fig^ The 
continuum production shown by the diagrams (a) and 
(b) in Fig'2] is often a background to direct resonance 
(p°, <fi, f2, etc.) production shown in diagram (c). We 
shall refer to (a)-(b) and (c) as the continuum and res- 
onance contributions, respectively. The zigzag line rep- 
resents the pomeron and subleading reggeon exchanges. 
The cross section for diffractive photoproduction of the 
opposite charge meson pairs, can be written as 
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da = (2tt) 6 (q + p — p+ — p~ — p 



d 3 p + d 3 p_ 



d 3 p' 
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flux 
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where for photoproduction, flux — A^J (p ■ q) 2 
2(s — m 2 ) and for a three-body reaction the amplitude 

[M] carries dimension of GeV - . The unpolarized cross 
section is a function of the following variables, the square 
of the center of mass energy, s = {q + p) 2 , the two-meson 
invariant mass, Mmm, ^mm = (^+ + ^-) 2 j the four- 
momentum transfer squared in the nucleon, t = (p — p 1 ) 2 
and the polar and azimuthal angles, £1 — (8, <fi) specifying 
the direction of momentum of one of the two produced 
mesons in the the two-meson recoil center of mass (RCM) 
system. The coordinates of the RCM system are usually 
chosen such that the y axis is perpendicular to the pro- 
duction plane determined by the photon and the recoil 
nucleon momenta, and z is chosen either in the direction 
of the photon momentum (in the RCM system), which 



is then referred to as the Gottfried- Jackson (GJ), or z 
is anti-parallel to the recoiling nucleon direction, which 
defines the so called s-channel helicity (SCH). At high 
energy (s M 2 IM ,t) the cross section takes a simple 
form 



da(M 



MM > 1 



(3 1 



dM MM dtdQ 



\ Ml p^M+M-p\2 j ( 2 ) 



where (3 



1 



( 2lL 



is the magnitude of velocity 



y Mmm ( 

of each meson in the center of mass of the M + M~ sys- 
tem. The invariant amplitude for the 2 — > 3 continuum 
process can be written in the Regge factorized form cor- 
responding to the diagrams (a) and (b) in Fig. 2] 



J 



V 7 m+ (A 7 ) 



Foe(t+) 



2 M^ p ( S -,t) + V lM -(^) 

l M 



F os (t-) .,M+ 



M" p (s + ,t). (3) 



'M 



In the equation above s+ and s_ are the Mandelstam 
variables for the M + p and M~p elastic scattering and 
t + and f_ are squares of momenta of the virtual mesons 
(t± = (k± — q) 2 )- For pseudoscalar mesons, because 
qe±(A 7 = ±1) = 0, the corresponding vertex functions 



read 



V lM ±{\ 1 ) = ±e (2fc£)e„(A 7 ) 



(4) 
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FIG. 1: Difltactive photoproduction of opposite-charge me- 
son pairs. The wavy line corresponds to the photon, the 
zigzag line describes the pomeron and subleading reggeon ex- 
changes, the lower lines are incoming and outgoing nucleons. 



subprocesses. At the high energies the invariant ampli- 
tudes for the 2^2 quasi-elastic subprocesses can be 
written in the simple but rather accurate form, 

M™>(s ± ,t) = is ± <ffi ±p (s ± ) exp (|t) 5\\> . (7) 

The Kronecker 5w reflects explicit imposition of the 
target nucleon helicity conservation, known to hold at 
high energies. The total cross section for irp and Kp 
are well known j7j and at high energies one can use the 
Donnachie-Landshoff parametrizations [2{|. For heavy 
mesons an educated guess will be necessary. In general, 
the heavier the mesons, the smaller the corresponding 
total cross section. The factor F os (t±) in Eq.® takes 
into account the extended nature of the exchanged par- 
ticle. Since for the process considered there are two ver- 
tices with an off-shell pseudoscalar meson, it is natural 
to write the combined form factor in the factorized form: 



F os (t±) = Fj% s (Q*,t±,m M ) ■ F corr (t ± ) 



The matrix elements take a simple form in the Gottfried- 
Jackson frame, 



e(A 7 = ±1) • k + = =F-^= sin(6> G j) exp(±i(f> G j) 
k 

e(A 7 = ±1) • fc_ = ± — sm(9 GJ ) exp(±i(f> G j) 
V2 



(5) 



The denominators in Eq.© can be calculated in terms 
of the canonical variables t, Mmm and cosd as 

t± - m 2 M = -2qk ± = ~(M MM - t)(l t PcosO) . (6) 

For small invariant masses Mm+m- of the two- meson 
state, (in the case of light meson production) sufficiently 
large overall s leads to large s+ and s_ in the relevant 



(8) 



where the first factor is the half-off-shell electromagnetic 
form factor from the upper vertex and the second one is a 
form factor representing the middle vertices in diagrams 
(a) and (b) in FigQ The exact form of the form factors 
is not known. In principle, a good quality data would 
help to find the proper functional form. 

Energy conservation imposes natural limits on energies 
in different two-body subsystems 



W + = ^/J^ <W — TTLM 

W- = < W - mm 
M M m <W -m p 



(9) 



The amplitude given in Eq.© is not yet complete since 
it does not satisfy electromagnetic current conservation. 
This current is given by, 



= 2c 



F os {t + ) M M- P{s _ t) 
+ t + -m\ [ 



Fos{t-) KA M+p 



t- 



(10) 



M 



so that the amplitude in Eq.® can be written as, 

e A1 (A 7 )J^(s,t, S+ ,t + ,s_,<_) (11) 



j^-yp^M+M-p 



Current conservation implies q^J^ = while from 



-4e 



F os (t + )M^ p (s_,t) 



-F os (t_)M™, p (s + ,t) 



(12) 



The origin of current non-conservation is two-fold. It 
comes form the non-point-like nature of the exchanged 



particles which introduces the from factors, F os ^ 1 
and from the difference in meson-nucleon scattering for 
the two charged mesons. The later implies that electro- 
magnetic charge flows differently in M + p — > M + p and 
M ~p — * M~p subprocesses and since photon couples to 
all charge currents there has to be a correction which 
reflects this difference. 

Having identified the two sources which contribute to 
the current we can unambiguously find the required cor- 
rections. We want to separate the corrections to the cur- 
rent which arise from interactions in the upper (meson) 
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and lower (baryon) vertices and therefore we define. 
1 



F ± = -[F os (t + )±F os (t-)] 



and 



M+ == - 



M + p, 



The current can be written as 



J" = 



7 M 

N 



where 



J c = 2e 



M 



(13) 



(14) 



(15) 



is a conserved current, 



(F+M+ + F-M-) 

(16) 



J^ = 2e 
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is non-conserved due to a difference between M + N and 
M~N cross-sections, and 



7 P — 2r 



t- 



t- 



''M 



(18) 



is non-conserved due the extended structure of the ex- 
changed meson. Now the additional contribution to the 
current required by current conservation will depend on 
meson variables for and nucleon variables for JjJ- , re- 
spectively, 



N 



N 



with 



and 



with 



2e 



q(p + p r ) 



F + M_ 



TP- 



TP- 



(19) 
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FIG. 2: Angular distribution of it in the Gottfried- Jackson 
frame. The dashed line is the contribution without corrections 
restoring current conservation. The dotted line represent the 
cross section associated with the corrections itself and the 
solid line corresponds to the full, current conserving ampli- 
tude. In this calculation we used the monopole off-shell form 
factor with M os = 1 GeV (see Eg. 152th 



is easy to check that the correction is negligible for the 
invariant mass distribution. 

In this paper we shall concentrate on the distribu- 
tions in Mmm and cos#. Often one is interested in 
distributions in the s-channel hclicity frame and not in 
the Gottfried-Jackson frame. From the dchnition of 
the two frames it follows that the spherical angles in 
the GJ frame can be expressed by those in the SCH 
frame M l P ^ M+M ~ p (t, M MM ,6 G j,<f> GJ ), with 9 GJ = 
" '"sch, 4>sch) and <j) GJ = <Pcj(0 S ch, 0scffl_through 



a rotation around the y axis by an angle 9 ro t 12 6l which 
in the high energy limit is given by, 



6 = 2e 



q(k+ + k-) 



F-M+ = 



^k ± ±k r Y F _ M+ 



t 



M 2 

MM 



(22) 

The the unphysical pole at t = M^ M , Ea.l|2*2^l should be 
eliminated by a zero in the F— form factor vanishes at 

t = M 2 MM- 

In Fig|2we show the effect of these corrections on the 
angular distribution for tt + it~ production calculated in 
the GJ frame. The correction is generally very small, 
except for the tips of the angular distributions. 

In this example we have integrated over the invariant 
mass range M T7r = 0.45-0.95 GeV and over the kinemat- 
ically accessible range of the momentum transfer t. It 



cos6> rot = 



M llM + t 



M MM 



t 



(23) 



In FigOl we compare angular distributions for the con- 
tinuum in the Gottfried- Jackson (left panel) and in the 
s-channel helicity (right panel) frames. The shapes in 
both frames are rather different. The effect of the rota- 
tion was neglected in the early calculations. 

In order for the Regge parametrization of the two-body 
amplitudes to be reliable, for 7r + 7r~ production the ener- 
gies W+ , W- have to be at least larger than 2 GeV. This 
can be fulfilled for center of mass energy W > 4 GeV, i.e. 
Ej > 8 GeV. The future GlueX experiment at TJNAF 
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is therefore at the border of application of the present 
model. 



III. RESULTS 



The cross section for a three-body reaction depends 
on five independent kinematical variables. For the reac- 
tions considered it is customary to use Mmm, and 
<j) and calculate <1(j(Mmm, t, fl)/dMMAidtdfl. The invari- 
ant mass distribution is then obtained by integrating over 
the remaining variables 
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da 



dM 



MM 



' M \un^ M ™l^\ (24) 



t m m{M M u) 



dMMMdtdVl 



where t m in and t max are calculated from the three-body 
kinematics (see for instance |3(J)- 



Light meson pairs 



FIG. 3: Angular distributions of ix + from continuum ampli- 
tude in the Gottfried- Jackson (left panel) and in the s-channel 
helicity (right panel) frames for four different energies: W = 
5 GeV (solid), W = 10 GeV (dashed), W = 20 GeV (dotted) 
and W = 70 GeV (dash-dotted). In this calculation we used 
the monopole off-shell form factor with M oa — 1 GeV (see 
Eq.JH2J). 



Below = 1 GeV the p° meson dominates the two- 
pion invariant mass spectrum. The amplitude for the 
resonance is taken in the relativistic Brcit-Wigner form. 
In our simple approach the normalization constant at t = 
is fixed based on the vector meson dominance model. 
We write the resonant three-body amplitude as, 



\s,t,M^,e,cj>) = C c 



7p 

r 



M p x £(s,t) f B w{M^)Yi,xJ6,< 



(25) 



Here we have introduced the amplitude for quasi-elastic 
scattering of p° meson off the proton 



iscjpopis) exp 



Bppt 



(26) 



As for the continuum model the Kronecker 5\\> reflects 
high-energy helicity conservation in the proton vertex. 
The factor C conv is adjusted to reproduce the correct 
normalization of the amplitude. 
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Different normalization of the reaction amplitudes are 
used in the literature. Throughout the present paper we 
use a popular in high-energy diffraction (see e.g. |24[) 
normalization, such that the angular distribution for the 
two-body reaction is 



da _ 1 / Vj_ \ jtj. — 12 

dn CM ~ 6A^ S U i lJVlfi] 



and the optical theorem at high energy reads 
ImM(s,t = 0) = sa tot {s) . 



(27) 



(28) 



This fixes the C CO nv factor in Ea. l|25JI . 

The factor fsw is the standard relativistic Brcit- 
Wigncr propagator, 



fBwiMirir) 



y/M T(M v „)/TT 

Ml 



- iM T(M„ x ) 



(29) 



where T(M 3r7r ) = T a ( 



mi 



-4m: 



M,?— 4m?. 



3/2 



and it is normalized 



according to $ \f B w{M^)\ 2 dM^ = 1 if T(M^j) is 
replaced by the energy independent width, IV We take 
M = 768 MeV and T = 151 MeV and 



tot 

a p°p 



(s) 



1 /> 

2 V" 



to L{s) 



,(«) 



(30) 



As for the continuum model the total cross sections for 
■n + p and n~p are taken from the Donnachie-Landshoff 
parametrization [25 ] . In the present approach we assume 
the same slope parameter for the resonance and contin- 



uum contributions, B a 



B. Except for the 



off-shell dependence determined by the form factors F os 
our model is essentially parameter free. We used different 
parametrizations the form factor: the exponential form, 



exp 



t± — m 



2AL 



and the monopole form 



F os {t±) 



Ml 



'M 



Ml - t± 



(31) 



(32) 



The form factor is normalized to unity at the on-shell 
point t = m 2 M . The exponential form is useful as the uni- 
versal parameter A os for different meson exchanges can 
be used. In the monopole parametrization M os > mm, 
i.e. different cut-off parameters for different exchanges 
have to be used. On the other hand the monopole form 
seems to be preferred, because at small virtualities vector 
dominance applies and in addition it produces a correct 
pQCD dependence at large virtualities. 

The off-shell form factor is the least know element of 
our model. In Figs0]|S]we present spectral shape (reso- 
nance + continuum) for the two choices of from factors 
without and with rotation from SCH to GJ frame. Some- 
what better agreement is obtained without performing 
the extra rotation of arguments. 
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FIG. 4: The spectrum of invariant mass M nlT for W = 70 GeV 
for exponential and monopole off-shell form factor without 
rotation of continuum arguments. The experimental results 
of the ZEUS collaboration are from 0. In this calculation B 
= 8 GeV" 2 . 



The spectral shape depends on the value of the form 
factor parameter (A os = 0.5, 1.0, 2.0 GeV for the expo- 
nential form factor, and M os — 0.5, 1.0, 2.0 GeV for the 
monopole form factor) . For pion production very similar 
results are obtained with exponential and monopole form 
factors for A os w M os . Therefore, having in view a pos- 
sible universality of the exponential off-shell form factor 
parameter A os , we shall use the exponential form factor 
in the following. The experimental data from the ZEUS 
collaboration at DESY [8j are superimposed on the the- 
oretical lines. In the absence of other mechanisms the 
value of the off-shell form factor could be obtained from 



7 



.100 





100 




90 




C 




\ 


80 




70 




60 


I 


50 




"O 
\ 


40 


b 


30 


u 






20 




10 








1 1 I 1 1 1 I 1 1 1 I 
W = 70 GeV 

exponential f.f 




.5 0.6 0.7 0.8 0.9 1 

M™ (GeV) 

(a) 



> 

o 



=1 



b 
u 



00 
90 
80 
70 
60 
50 
40 
30 
20 
10 




1 1 i 1 1 1 i 1 1 1 i 
W = 70 GeV 

monopole f.f 




0.5 0.6 0.7 0.8 0.9 1 

NC (GeV) 

(b) 

FIG. 5: Same as in Fig. [1] but with rotation of continuum 
arguments. 
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FIG. 6: The spectrum of invariant mass M^^ for W = 70 
GeV. The experimental results of the ZEUS collaboration arc 
from In this calculation B = 8 GeV -2 . The standard 
resonance contribution is shown as a reference (dashed line). 



effect which deforms the spectral shape of the p meson. 

In our approach the background is described by a 
physical process which has a strong two-pion invari- 
ant mass dependence. Often in experimental analyses 
IE El US (see also this background is parametrized 

with a weak irir mass dependence. We do not expect 
that our model could describe the experimental spec- 
tra above M W1T > 1 GeV, since there will be important 
contributions from higher mass two-pion resonances e.g. 
/ 2 (1270), p 3 (1690), ,9(1700). The situation there may 
be therefore rather complicated and we leave the corre- 
sponding analysis for future investigations. 

In order to further understand the large interference 
effect of the continuum contribution with the P-wave res- 
onance we performed a decomposition of the continuum 
amplitude into the partial wave series in the GJ frame, 



the fit to the experimental data. The coherent sum of 
the resonance and continuum (solid line) differs consid- 
erably from the standard resonance shape which is shown 
in FigEl The resonance contribution alone (dashed line) 
gives a poor description of the data. In particular, the 
position of the maximum is at higher m v „ then observed 
experimentally. These features have been often ignored 
in the literature and only integrated cross section were 
used to compare with theoretical calculations. 

In order to have a better insight into the origins of the 
line-shape modifications in Fig[7|we show separately the 
resonance and the continuum contributions in somewhat 
broader range of two-pion invariant mass. This figure 
clearly demonstrates that it is mainly the interference 



Ml^ + x r p {t,M^6,<t>) = 

Lin 

(33) 

The expansion coefficients can be calculated as 

a hn ' (M'^tt) = 

= J Y? m (e, 0) • MY- x A S x M ~ p (t, M^,6, 0) rfO. 

(34) 
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FIG. 7: The decomposition of the spectrum of the ir + -ir~ 
invariant mass M W7I - for W = 70 GeV. We show separately 
the resonance (dashed line) and continuum (dotted line) cross 
sections. The solid line corresponds to a coherent sum of both 
processes processes. In this calculation B = 8 GeV -2 . 



In our model the expansion coefficients depend only on 
A 7 , 1, m, i.e. 



{t,M^) = aZ(t,M^) . 



(35) 



We find that in the case of the pseudoscalar production 
the continuum contributes dominantly to the P-wave, i.e. 
I aim | > |aoo|, |a 2m |, etc. and \af\\ < \a^\ < (a^l.and 

l a l~-il > \ a ia\ > l a i+il- This explains the large interfer- 
ence between the continuum and resonance production. 
We also find a relatively large contribution of the -F-wave. 
The individual contributions of I = 1 and I = 3 partial 
wave are shown in FigEl 

Presence of the F-wave is interesting in the context of 
forward-backward asymmetry and the moment analysis 
(see e.g.0). These are usually discussed in terms of the 
S- and P- wave interferences. Our analysis shows that in 
principle one needs to include S, P, D and F waves into 
such an analysis even for relatively low invariant masses. 

Let us try to understand this hierarchy of the par- 
tial wave amplitudes. The angular distribution originates 
from 
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72>— >M + M p 



(t, M M m, e,<f>) 

oc qMMi^(t) exp(±i<j)) sm6A(9, < 



FIG. 8: The contributions of individual partial waves to the 
spectrum of invariant mass M^-k for W = 70 GeV (excluding 
resonance production). In this calculation B = 8 GeV -2 . 



where we have defined a slowly changing function of 9 
and <t>. 



A{0,< 



a M - P (s-(e^))F os (t + (8)) 

1-/3 cos e 
a M+p (s + (O,<l>))F oa (t-(0)) 

1 + P cos 6 



(37) 



It then becomes obvious that it is the function A(6, <j>) 
which is responsible for generation of partial waves differ- 
ent than 1 = 1. It is easy to show that if the numerators 
were identical the S'-wave amplitude would vanish. The 
smooth energy dependence of the cross sections intro- 
duces a small S'-wave contribution. This small effect is 
slightly dependent on the incident energy. 

The large interference effect is specific to the photopro- 
duction of P-wave resonances and the P-wave dominated 
pseudoscalar meson continuum. 

In Fig|5| we discuss the evolution of the spectral shape 
(asymmetry) as a function of the momentum transfer, \t\. 
While at low \t\ the spectral asymmetry is reversed com- 
pared to the standard resonant shape at large \t\ one ob- 
serves a restoration of the standard asymmetry expected 
for the P-wave resonance. Such an effect was observed 
experimentally in [T(il | and to the best of our knowledge 
the dynamics of this effect has not been given before. 

So far we have neglected the final state 7T7t interaction. 
This can be restored by modifying the continuum partial- 
wave amplitudes, 



= CL,1' ' 



(36) 



:t XA ' cos SreM^D- (38) 

Here 5™ is the M W7r -dependent phase shift for 7T7t scat- 
tering in the Vth partial wave. If only resonant final state 
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FIG. 9: The spectrum of invariant mass M nn for W = 70 
GeV for different bins in t. The solid lines correspond to the 
results without and the dashed line to the results with the 
rotation of arguments of the continuum contribution. In this 
calculation B = 8 GeV~ 2 and exponential off-shell form factor 
with A = 1 GeV has been used. 



FIG. 10: The spectrum of invariant mass M K g. The ex- 
perimental results of the FOCUS collaboration (without ac- 
ceptance corrections) are from [3lll . The experimental data is 
not normalized. In this calculations B — 6 GeV~ 2 and and 
exponential off-shell form factor with A = 0.5, 1.0, 2.0 GeV 
has been used. 



interaction effects are included the modified partial am- 
plitude of the continuum becomes, 



~A-^,A,A A-^,A,A 
a lm = a lm 



Ml 



Ml 



M 2 



Ml 



iMqT 



(39) 



Such a modified amplitude vanishes at M v „ = Mq. For 
I = 1 the final amplitude can be written as a sum of three 
terms: the direct p° production, free meson pair pro- 
duction and the resonance production via re-scattering. 
The last two are contained in Eq.J^ (or Eq.lJSSJl). We 
checked that re-scattering from the continuum back to 
the resonance is negligible. 

The K + K~ production is more sensitive to the choice 
of the off-shell form factor. The spectrum of two-kaon 
invariant mass is shown in Fig 1 101 for exponential form 
factor and A os = 0.5, 1.0, 2.0 GeV. 

Recent experimental data from the FOCUS collabora- 
tion are shown for comparison |l5j | . Our model of diffrac- 
tive production of kaonic pairs gives a good description of 
the main trend of the data. A direct comparison with the 
FOCUS data is, however, not possible, because the exper- 
imental acceptance is known only in very limited range of 
the phase space. The absolutely normalized experimen- 
tal invariant mass distribution would help to better limit 
the only unknown parameter A os (or M os ) of the off-shell 
form factor. The situation may be, however, not so sim- 
ple in light of resent experimental data for 77 collisions 
|l7|. In principle, similarly as for the 7r + 7r~ invariant 
mass the interference of the K + K~ continuum and the 
(j> resonance takes place. However, because the <p has a 



much smaller decay width than the p° and is situated 
very close to the K + K~ threshold (where the contin- 
uum contribution is small), its importance for Mkk > 
1.1 GeV is negligible. 

The angular distributions depend on the invariant 
mass. Thus one studies angular distributions in bins of 
invariant masses. The angular distribution is calculated 
from the 4-dimensional differential cross section as, 



da 



dcosO 



dt 



da(M M M,t,n) 



. dMMM j t . -j-- dM MM dtdn 

(40) 

In Fig^Jwe present angular distribution for the ZEUS 
kinematics 0. The angular distribution is almost pro- 
portional to (sin0) . The continuum contribution only 
slightly modifies the (sin#) 2 resonance distribution. The 
results with or without rotation between the frames are 
almost indistinguishable. 

In Fig ll2l we present predictions for angular distribu- 
tions for the GlueX experiment at TJNAF. The shape of 
the distributions depends on the interval of the two-pion 
invariant mass. The effect of the rotation of the argu- 
ments of the continuum is stronger for smaller invariant 
masses. We predict a sizable asymmetry with, n + being 
preferentially emitted in the forward (photon) direction. 

The function A(9,<p) in Ea. (|36() is responsible not 
only for generating higher partial waves but also for the 
forward-backward asymmetry. In order to measure the 
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FIG. 11: The angular distribution of charge pions in the 
SCH frame for W = 70 GeV and M-mr in the measured range. 
The experimental results of the ZEUS collaboration are from 
0. In this calculation B = 8 GeV -2 and exponential off-shell 
form factor with A = 1 GeV have been used. The solid line 
corresponds to the calculation without rotation of arguments 
of the continuum, whereas the dashed line corresponds to the 
calculation with the rotation. 



asymmetry we define the following quantity, 



By construction 



(41) 



(42) 



which means that the asymmetry must fulfil the symme- 
try relations 



A* FB (6) = -Af B (6) 



(43) 



The asymmetry for tt + is shown in Fig ll3l for four differ- 
ent energies and for one bin in M„, 0.55 GeV < M n7r < 
0.95 GeV and \t\ < 1 GeV 2 . We present separately the 
asymmetry of the continuum contribution alone (panel 
a) and the asymmetry of the sum of the resonance and 
continuum contributions (panel b). A sizable asymme- 
try can be seen in the GJ frame (solid lines). In general, 
the larger incident energy, the smaller the asymmetry. 
The inclusion of the resonance contribution lowers the 
asymmetry around z = 0. 

For the GlueX experiment at TJNAF we predict the 
asymmetry to be about 10% at the very forward and very 



FIG. 12: The angular distribution of ir + in the SCH frame 
for the GlueX energy W = 4 GeV. In this calculation B = 
6 GeV -2 and exponential off-shell form factor with A = 1 
GeV have been used. The sold line corresponds to the case of 
no extra rotation while dashed line corresponds to the result 
with the extra rotation. 



backward directions. Even at W = 70 GeV the asymme- 
try is of the order of 1%. In our model the asymme- 
try is caused by the different interaction of ir + p and ir~p 
(analogously for K + p and K~p). This is caused by a dif- 
ferent strength of sub-leading reggeons. The asymmetry 
disappears at large energy where the dynamics of elas- 
tic scattering is governed exclusively by the pomeron ex- 
change. The asymmetry discussed above is an important 
test of the diffractive mechanism. When transformed to 
the SCH frame the asymmetry becomes rather negligible 
(dashed lines). 

In Refs.[33, H^] the observation of forward-backward 
asymmetry of charged pions was proposed in order to pin 
down the odderon exchange. These analyses were based 
on the assumption that only resonant mechanisms plays 
a role. The mechanism considered here was not taken 
into account. Our diffractive mechanism may mimic 
the pomeron-odderon interference effects discussed in 
Refs.|32l l33|. A careful search for the odderon exchange 
must therefore necessarily include the two-pion contin- 
uum discussed in the present paper. 

We expect that the diffractive production of K + K~ is 
the dominant mechanism well above the <fi resonance. In 
FigEl we present angular distributions of K + in the GJ 
recoil center of mass system for a typical FOCUS energy, 
W = 10 GeV for different bins of Mkk specified in the 
figure. 

As for the pion production we obtain asymmetric dis- 
tributions with K + being preferentially emitted in the 
forward hemisphere. The shape of the angular distribu- 
tion changes with Mkk- The shape is very important 
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FIG. 13: Forward-backward asymmetry for n for the con- 
tinuum (panel a) and for the resonance+continuum (panel b) 
as the function of cosf? for four different energies: W = 4, 10, 
20, 70 GeV. The solid line corresponds to the continuum in 
the GJ, while the dashed line to the continuum in the SCH 
frame. The kinematical cuts are specified in the text. 



FIG. 14: The angular distribution of K in the GJ frame for 
a typical FOCUS energy W = 10 GeV. In this calculation B 
= 6 GeV" 2 . 




when studying resonances. For example, the FOCUS 
collaboration found difficulties in spin assignment of the 
X(1750) 1^3 . This may partially be due to the inter- 
ference of the resonance and the continuum. In general, 
the larger the invariant mass, the more the cross section 
peaks at forward/backward directions. This is related to 
the increase in the number of active partial waves with 
the increasing subsystem energy. 

In Fig E| we present Afb for K + for four different 
incident energies, 1.9 GeV < M KK < 2.1 GeV and |t| < 1 
GeV 2 . We observe much larger asymmetries than for the 
7r + 7r~ case. This is partially due to larger asymmetries 
in K + p and K~p scattering than for the ir + p and 7r~ p 
scattering. 

In principle the K + K~ asymmetry should be seen 
in the data of the FOCUS collaboration. The possi- 
ble higher mass resonance contributions are expected to 
lower the asymmetry. A careful search for the asymme- 
try as a function of the two-kaon invariant mass would 
be very interesting in searches for new states. 



FIG. 15: Forward-backward asymmetry for K + production in 
the GJ frame as the function of cos# for four different energies: 
W = 4, 10, 20, 70 GeV. The kinematical cuts are specified in 
the text. 



B. Heavy meson pairs 

In Fie llbl we show invariant mass distributions of 
diffractively produced pairs of D + D~ for the average en- 
ergy of the FOCUS experiment 22]. 

The averaging is not so important as the diffractive 
contribiution is only slowly dependent on the incident 
energy. We show results for different values of the scale 
parameter in the exponential off-shell form factor. The 
absolute normalization strongly depends on the value of 
the form factor mass parameter. This is due to the fact 
that the intermediate meson D in Fig^ is far from its 
mass shell. Our distributions peak at Mud = 4.5-5 
GeV. In principle, the FOCUS collaboration could an- 
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FIG. 16: The spectrum of invariant mass Mod- The theoret- 
ical results were obtained with B = 6 GeV -2 and exponential 
off-shell form factor with parameters specified in the figure. 



FIG. 17: The spectrum of invariant mass Mbb- The theoret- 
ical results were obtained with B = 6 GeV -2 and exponential 
off-shell form factor with parameters specified in the figure. 



alyze their data and try to construct the invariant mass 
distribution. The expected statistics will be of course 
rather low, of the order of 10-50 events. 

The invariant mass distribution for the B + B~ pair 
production is shown in FigEI f° r a typical HERA en- 
ergy W = 100 GeV. Here there is a dramatic effect on 
the value of the form factor parameter. The distribution 
reaches its maximum at rather high Mbb- The small- 
ness of the cross sections precludes, however, studies of 
differential cross section. 

In FigsElEHwe show the energy dependence of the 
integrated cross section for both DD and BB pair pro- 
duction (thick solid line). For comparison we show the 
standard collinear-factorization results for open charm 
and open bottom production [34[ as well as the experi- 
mental cross sections. As can be seen from the figure the 
diffractive cross section is much smaller than its coun- 
terpart for the standard photon-gluon fusion followed by 
fragmentation. 

In addition, the rise of the cross section with energy 
is slower then for the photon-gluon fusion. Only close 
to the kinematical threshold the diffractive component 
may constitute a sizable fraction of the open heavy-flavor 
component. Since in this case application of the meson- 
exchange approach is less reliable than for the light pairs, 
an experimental measurement would be very interesting 
and helpful to discriminate between models. While for 
the D + D~~ pairs this may be feasible, it will be hard 
to expect such an analysis for the B + B~ pairs, at least 
in a few year perspective. We predict the fraction of 
diffractive-to-nondiffractive events to be somewhat larger 
for bottom than for charm mesons. 



y + "p — > cT cbor 




1 10 1 o 2 

W (GeV) 



FIG. 18: The energy dependence of the integrated cross sec- 
tion for the diffractive photoproduction of D + D~ (thick solid 
line) with our estimate of the upper limit (A os = 2 GeV). 
For comparison we show the cross section for open charm 
and bottom production (dash-dotted line) and corresponding 
experimental data. The details concer ning the open heavy- 
flavour production can be found in Ref. |34f . 



IV. CONCLUSIONS 

We have presented a model of diffractive photoproduc- 
tion of opposite-charge pseudoscalar mesons. The model 
should be valid at sufficiently high energies, down to en- 
ergies relevant at the future experiments in Hall D at 
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FIG. 19: Same as in Fig.EHlfor B + B~ production, (A os = 
4 GeV). 



Jefferson Lab. The understanding of the two-meson con- 
tinua is absolutely crucial when looking for new (exotic) 
mesonic states in the two-meson channels. Our model re- 
sults seem to be consistent with the current experimental 
data. 

The interference of the two-pion continuum and the 
resonant p a contribution leads to a significant deforma- 
tion of the resonance peak, which is consistent with 
the experimental data from the ZEUS collaboration at 
HERA 8]. The effect of the continuum-resonance in- 
terference is often completely neglected. We explain a 
change of the spectral shape of the p° bump with mo- 
mentum transfer. 

The diffractive mechanism of the K + K~ produc- 
tion leads to a broad bump in invariant mass distribu- 
tion, above the 4> meson position with the maximum at 



M K k ~ 1.4 GeV. Experimentally, except for the X(1750) 
mentioned earlier no clear resonances in the K + K~ chan- 
nel have been seen in |l5| e.g. even though states like 02, 
fi, P3 are expected. A careful phase-shift analysis will 
be required to separate them. Our model amplitude pro- 
vides a reasonably well controlled background for such 
studies. 

We have found a forward-backward asymmetry in n + 
and 7r~ (similarly K + and K~) recoil center of mass 
angular distributions. The new effect is due to differ- 
ent interaction of different-charge mesons with the pro- 
ton. The effect disappears at large incident energies, 
where the pomeron exchange is dominant. This asym- 
metry may be very important in the context of recently 
proposed searches for odderon exchange via forward- 
backward charge asymmetry in the region of M nv <~ 1 
- 2 GeV. 

The cross sections for diffractive production of heavy 
meson pairs were neither calculated nor measured in the 
past. In the present paper we have presented an esti- 
mates of the corresponding cross sections. The cross 
sections obtained in our analysis are, however, very 
small. Even a measurement of integrated cross section for 
diffractive D + D~ pair photoproduction would be very 
useful for testing the mechanism of diffractive produc- 
tion. Such data would allow one to pin down off-shell 
effects for the intermediate (exchanged) heavy meson. 



Acknowledgments 

We thank Ryan Mitchell for providing us the experi- 
mental two-kaon invariant mass distribution of the FO- 
CUS collaboration. We are also indebted to Leonard 
Lesniak and Geoffrey Fox for several discussions. This 
work was partially supported by the US Department of 
Energy grant DE-FG0287ER40365 (APS), and the Nu- 
clear Theory Center at Indiana University. 



[1] S.D. Drell, Phys. Rev. Lett. 5, 278 (1960). 

[2] P. Soding, Phys. Lett. 19, 702 (1965). 

[3] A.S. Krass, Phys. Rev. 159, 1496 (1967). 

[4] G. Kramer and J.L. Uretsky, Phys. Rev. 181, 181 (1969). 

[5] J. Pumplin, Phys. Rev. D2, 1859 (1970) . 

[6] G. Kramer and HR. Quinn, Nucl. Phys. B27, 77 (1971). 

[7] Review of Particle Physics, Particle Data Group, Eur. 

Phys. J. C15, 1 (2000). 
[8] M. Derrick et al. (ZEUS collaboration), Z. Phys. C69, 39 

(1995). 

[9] J. Breitweg et al. (ZEUS collaboration), Eur. Phys. J. 

C6, 603 (1999). 
[10] J. Breitweg et al. (ZEUS collaboration), Eur. Phys. J. 

C14, 213 (2000). 
[11] DR. Thompson et al. (E852 collaboration), Phys. Rev. 

Lett. 79, 1630 (1997). 



[12] A. Abele et al. (Crystal Barrel collaboration), Phys. Lett. 
B423, 175 (1998). 

[13] A. R. Dzierba et al, Phys. Rev. D 67, 094015 (2003). 

[14] A. P. Szczepaniak, M. Swat, A. R. Dzierba and S. Teige, 
Phys. Rev. Lett. 91, 092002 (2003). 

[15] J.M. Link et al. (FOCUS collaboration), Phys. Lett. 
B545, 50 (2002). 

[16] A. Sandacz (COMPASS collaboration), private commu- 
nication. 

[17] K. Abe et al. (BELLE collaboration), Eur. Phys. J. C32, 
323 (2003). 

[18] Ch-R. Ji, R. Kamihski, L. Lesniak, A. Szczepaniak and 
R. Williams, Phys. Rev. C58, 1205 (1998). 

[19] L. Bibrzycki, L. Lesniak and A. P. Szczepaniak, Eur. 
Phys. J. C34, 335 (2004). 

[20] S. Chekanov et al. (ZEUS collaboration), 



14 



hep-ex/0405069 
[21] A. Csilling (OPAL collaboration), in PHOTON2004, 
Ambleside, UK, August 26-31, 2000, hep-ex/0011022 
M. Acciari et al. (L3 collaboration), Phys. Lett. B503, 
10 (2001); 

F. Kapusta (DELPHI collaboration), in MESON2004, 

Cracow, Poland, June 4-8, 2004. 
[22] J.M. Link et al. (FOCUS collaboration), Phys. Lett. 

B566, 51 (2003). 
[23] M. Luszczak and A. Szczurek, Phys. Lett. B594, 291 

(2004). 

[24] V. Barone and E. Predazzi, it High-Energy Particle 

Diffraction, (Springer, Berlin 2002). 
[25] A. Donnachie and P.V. Landshoff, Phys. Lett. B296, 227 

(1992). 

[26] K. Schilling, P. Seyboth and G. E. Wolf, Nucl. Phys. B 



15, 397 (1970) [Erratum-ibid. B 18, 332 (1970)]. 
[27] M.G. Ryskin and Yu.M. Shabelski, Phys. Atom. Nucl. 

61, 81 (1998). 
[28] D.C. Fries et al., Nucl. Phys. B143, 408 (1978). 
[29] D.P. Barber et al., Z. Phys. D12, 1 (1982). 
[30] E. Byckling and K. Kajantie, Particle Kinematics, (John 

Wiley and Sons, London 1973). 
[31] R.E. Mitchell, Phd thesis, the University of Tennessee, 

Knoxville, December 2003. 
[32] Ph. Hagler, B. Pire, L. Szymanowski and O.V. Teryaev, 

Eur. Phys. J. C26, 261 (2002). 
[33] I.F. Ginzburg, LP. Ivanov and N.N. Nikolaev, Eur. Phys. 

J. C30, 2 (2003). 
[34] A. Szczurek, Eur. Phys. J. C26, 183 (2003). 



